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Abstract 

In this paper we consider the strong asymptotic behavior of Laguerre polynomials in the complex 
plane. The leading behavior is well known from Perron and Mchlcr-Hcine formulas, but higher 
order coefficients, which are important in the context of Krall-Laguerrc or Lagucrrc-Sobolev-typc 
orthogonal polynomials, are notoriously difficult to compute. In this paper, we propose the use of an 
alternative expansion, due to Buchholz, in terms of Bessel functions of the first kind. The coefficients 
in this expansion can be obtained in a straightforward way using symbolic computation. As an 
application, we derive extra terms in the asymptotic expansion of ratios of Laguerre polynomials 
in C\ [0,oo). 

AMS subject classification (2010): 33C45, 30E15, 33C10 

Keywords and phrases: Laguerre orthogonal polynomials, asymptotic expansions. 

1 Introduction 

The classical Laguerre polynomials (x) are defined as the polynomials orthogonal with respect to 
the 1? inner product 

oo 



(P,q) a = / p(x)q(x)x a e x dx, a > -1, p,q€ 
Jo 

see for instance [6] or [16]. We normalize them according to j!6j : 



(_1)n 

L£Hx)= { —LLtHx), (l) 
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where iffl (x) denotes the monic Laguerre polynomial of degree n. 

These polynomials can be given in terms of an \F\ confluent hypergeometric function as follows 
(see, for instance, p], [16]): 

L£\x)= ( n + a ) 1 F 1 (-n;a + l;x). (2) 

The asymptotic behavior of Laguerre polynomials in the complex plane as n — > oo is well known 
in the literature. In particular, we recall the following two asymptotic expansions: 

• Outer strong asymptotics (in C \ ]R + ) is given by Perron's formula. For a > — 1 we get 

4") ( z ) = -}—el 2 (-^)-«/ 2 -V4 n a/2-i/4 e2v ^ | £ Cm (z)n- m / 2 + 0(n~ d / 2 ) } . (3) 



,m=0 



Here Ck{z) is independent of n. This relation holds for 2 in the complex plane with a cut along 
the positive real semiaxis. The bound for the remainder holds uniformly in every closed domain 
of the complex plane with empty intersection with IR + (see |16j . Theorem 8.22.3). 

• Mehler-Heine type formula. Fixed j £ Z + and J a the Bessel function of the first kind, then 

lim L ^M^±M = z -^ Ja ( 2 ^) , (4) 

uniformly over compact subsets of C (see |16j . Theorem 8.1.3). 

Higher order coefficients in the asymptotic expansion are important when one deals with Krall- 
Laguerre or Laguerre-Sobolev-type orthogonal polynomials. More precisely, they play a key role in 
the analysis of their outer relative asymptotics, see [9]. There one needs to estimate ratios of Laguerre 
orthogonal polynomials like 

where n = 0, 1,2, . . ., j £ Z. Additionally, we require a, /3 > —1. Note that in the particular case 
where j is an integer and a = f3, this asymptotic information can in principle be obtained from the 
three-term recurrence relation for Laguerre polynomials, applying the Perron theorem, see for example 
[TUl §4.3], but for general values of a, (3 and j this procedure is not feasible. 
We consider the Laguerre kernel 

„ , , n + 1 Li%(x)L < n a \y)-Lt\x)L^ l (y) 
K n -i{x,y) = — ^y- 



and we denote by 



in u J 1 12 x — y 



K^(x,y) = 



dxWy k 

their partial derivatives of order j and k with respect to the variables x and y, respectively. 
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If we compute the d-th. derivative with respect to the second variable and we evaluate it at y = c, 
that is Kn' d \x,c), we find expressions of the form 



d 

£ 

3=0 



, )i+ l L< n-(l+d)+j ( C ) 



(5) 



For instance, when d = 1 we get 



n\(n-l)\ r ( a ) , v ( a +i) 



llfW 112 

,4- + 2 1} ( C ) 



n-2 




^ a) (x) ' 



+ 



(6) 



see 0. When d = 2, we get 



i 



n-2! 

r(«) 112 
(a+2) 



(a+2) 
n— 3 



n-1 



(c) 



(x — c) 3 

1 



(a+2) 
n— 3 



(c) 



'4- + 2 2) (c) 



(x 



(a+2) 
n— 3 




(7) 



3 » 4 a A(*) 



These formulas follow from differentiation and using two well known identities for Laguerre poly- 
nomials. The first one is the so called structure relation 



LW(x) = 4° +1) (*) - (*) , n > 1. 

The second one is the Hahn property 

[L^(x)]' = -L^\x), n>l, 



(8) 



(9) 



see for example [16], formulas (5.1.13) and (5.1.14). 

When considering the asymptotic behavior for n large enough, one needs higher order coefficients 
in n for those terms of the form ([5]) in the previous expression, because of cancellation. More precisely, 
we need to know exactly the coefficient of n~ d ^ 2 to estimate ([5]) correctly. For example, if d = 1 we 
need to know the coefficient of n _1//2 , if d = 2 the coefficient of n _1 , and so on. Thus, expressions like 



r (o+i) 

lim „NAV*J 



L 



n+h 



(x) 



(-x) 



h,k S Z, 



appearing in [21 §8], or Lemmas 1 and 2 in [9], are not accurate enough if d > 3. 

The main advantage of Perron's expansion for Laguerre polynomials is the simplicity of the asymp- 
totic sequence (inverse powers of n), but it has the problem that the coefficients Ck(z) soon become 
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cumbersome to compute. The standard expansion, valid for z outside of the support of the orthog- 
onality measure, is due to Perron [T5], based on ideas of the steepest descent method, and then it 
appears referenced in [16] and many others. 

One possibility is to use the generating function for Laguerre polynomials: 



write the coefficients as contour integrals and apply the standard method of steepest descent, see for 
instance [3] or [12]. However, the computations soon become complicated, since parametrizing the 
path of steepest descent is not easy in explicit form. 

In principle, it is also possible to compute higher order terms in Perron's expansion using Riemann- 
Hilbert techniques. Following ideas of [7], in [20] the author presents strong asymptotics of orthogonal 
polynomials with respect to the weight w(x) = e~®^ on [0, oo), where Q(x) is a general polynomial 
with positive leading coefficient. This clearly includes the Laguerre case. In order to obtain higher 
order terms, it is necessary to work further into the last matrix of the nonlinear steepest descent 
method, R(z), and then point out the details undoing the transformations. 

To the best of our knowledge, the only sources of information for higher order coefficients in the 
Perron expansion are the works by Van Assche [181 119] . and D. Borwein, J. M. Borwein and R. E. 
Crandall in [3J. This last one is based on complex integral representations with strict error bounds. 
It provides a powerful method to generate the coefficients C m . 

In this paper we propose an alternative approach, based solely on using an expansion of the 
Laguerre polynomials that involves more complicated special functions, namely Bessel functions of 
the first kind. This type of expansions go back to the works of Tricomi and Buchholz, see [5j §7.4], 
and also [15] . [HJ §10.15] and references therein. In this way, the different behaviors of Ln{%) in the 
complex plane are better captured, and thus the coefficients are simpler. Moreover, apart from the large 
n asymptotic behavior, the resulting approximation converges in the complex plane. Reexpanding the 
Bessel functions in inverse powers of n destroys this property, but using straightforward manipulations 
it is possible to recover the results in [3]. 

The structure of the manuscript is as follows. In Section 2 we analyze some expansions of confluent 
hyper geometric functions in terms of Bessel functions of the first kind. Thus we deduce an algorithm 
which gives the polynomial coefficients of such expansions for Laguerre polynomials. This approach 
has many advantages with respect to the standard Perron expansion. Based on this method, in 
Section 3 we deduce the asymptotic expansions of outer ratio asymptotics for Laguerre polynomials 
which cover higher order terms than those used in the literature until now. 

2 Tricomi and Buchholz expansions in terms of Bessel functions 

In [IT], the authors study several expansions of confluent hyper geometric functions in terms of Bessel 
functions, introduced by Tricomi and Buchholz. Some of these expansions are useful in an asymptotic 
sense. We have 




z\ < 1 




(10) 



771=0 
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where 

K =--a (11) 

is the Whittaker parameter, the function E a {z) can be expressed in terms of the standard Bessel 
function of the first kind 

E a {z) = z- a / 2 J a (2^), (12) 
and the coefficients A n (a, c) satisfy the recursion 

(m + l)A m+1 = (m + c- l)^ m _i - 2kA ti _ 2 , m>2, 

with initial values 

c 1 
A (a,c) = l, Ai(a,c) = 0, A 2 (a,c) = -, A 3 (a, c) = --(c - 2a). 



The expansion (jlOp is convergent in the entire z plane. This follows from results of Tricomi, see 
IE and [m. 



For Laguerre polynomials L^\z), we set a = —n and c = a + 1, i.e. 

4 Q) (*)= M-ma + liz). 

Observe that for positive real values of z, the Bessel functions J a {2y/nz) are oscillatory, a fact 
that is consistent with the behavior of the Laguerre polynomial in the interval of orthogonality. Also, 
because of the construction of the function E a (z), the series expansion is also asymptotic for large 
values of k. Indeed, k = n + so the argument of the Bessel function becomes large when the 
degree of the Laguerre polynomial n grows. Note also that if we take the first term in this expansion 
and scale the variable z, we recover the Mehler-Heine asymptotics However, we point out an 
important drawback of this expansion: the coefficients A n {a,c) are polynomials in a, and hence they 
grow with n in the Laguerre case. 

An alternative expansion given in [11] is due to Buchholz, see §7.4], and has a similar structure: 

oo 

l F l {a-c-z)=T{c)e z l 2 ^ (g) P m (c,z)E c _ 1+m (Kz), (13) 

771=0 

where n is introduced in (fTTj) and E c _i +m (Kz) is given in (fl~2|) . The new coefficients P m (c,z) satisfy 
a more complicated relation: 

P m (c,z) = z~ m l 2 J* (^uP m ^(c,u) + {c-2)P' m _ x {c,u) - uP^^u)) vT^du (14) 

for m > 1, with initial value Po(c, z) = 1. The first coefficients are 

P 1 (c,z) = ^, P 2 (c,z) = ±-(z 2 + l2c-24), P 3 (c,z) = -A-(5z 2 + 180c - 432). 
b 72 6480 

An important advantage of this expansion is that the coefficients do no longer depend on a, so 
they remain bounded in the Laguerre case as n — > oo. Actually, as explained in [11], the expansion 
has an asymptotic character for large values of k (i.e. large values of n) with kz > 0, when \z\ = o(k). 
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Thus, we obtain 



L^(z) = r + a J r(a + l)e z ' 2 (| ) Pm(a + 1, z)E m+a (Kz) 

^ ' m=0 

r(n + a + 1) /2 ^ (Z\ m . w s_rri+a T , n , — s 

= — j -e ' Z^{o) p m(a + l,z)(Kz) 2 J m+Q (2^) 



n! 

m=0 



r(n + a + i) eZ/2M _ a/2 £ / M-/ 2 Pm(a + M)Jm+a(2 ^ ) . 



(15) 



n! ' \4m 

m=0 

Again it is possible to recover the Mehler-Heine asymptotics from this expression. 

The aim of this section is to explain how to construct an asymptotic expansion in terms of negative 
powers of n from (|15p in a systematic way. The coefficients become fairly complicated as well, but the 
procedure is easily implemented using symbolic computation. 

In order to rewrite (|15p in terms of negative powers of n, we use the asymptotic approximation for 
the Bessel function of large argument, [TJ 10.17.3]: 

V 7 V fc=0 k=0 / 

where _ 

u = u(a) = z- — - - (16) 

and the coefficients are 

, - (4a 2 - l)(4a 2 - 9) . . . (4a 2 - (2fc - l) 2 ) 
a k (a) = ^ L_A i L± t k > 1, (17) 

ao(a) = 1. Consequently, 

/no , oci 

k a 2 k+i[m + a) 



(00 / . 00 

cos^^(-l) fc a2fc ^ + Q) - S ina,]T(-l)' 
k=0 ^ KZ > k=0 



^ KZ )k+l/2 J ' 

for j arg(/-tz)| < 2tt and integer m > 0. Using these expressions, we obtain the following result: 

Theorem 1 For a > — 1, the Laguerre polynomial L^*\z) admits the following asymptotic expansion 
as n — » 00 : 



L() = r(n + a + l) 7r _ 1/2e , /2 /2 _ 1/4 



d „ d 



E B2m.(a,z) cos uj ^ B 2m +i{a,z)smu} 



m=0 m=0 



(18) 



/or some coefficients B m (a, z) independent of n. The error term holds uniformly for z in compact sets 
of C, and the parameter uj = 00(a) is given by (|16p . 
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Here, the coefficients B m (a,z) come from the expansion of the Bessel functions. We observe that 
we have a sum of Bessel functions of different orders, but the sines and cosines that appear in the 
asymptotic expansions can be somehow grouped together, since 



cos ( z - {a + m)n -l)=(-iy\ cos (* - " - i) ' m = 2s ' 



cos 


(*" 


Q7T 
2 


TT 
1 


sin 


(*- 


CK7T 


7T~ 


2 




sin 


(*" 


Q7T 


7T 


2 





sin [ z _ (" + m ^ -l) = (-iy J sin ( z - T - f ) ' m = 2s ' 

2 4y 1 -cos(-g- ^ - f) , m = 2s + l 



(19) 



for s = 0,1,2, ... 

An extra care is needed in two aspects when assembling the expansion, though: first, one has to 
take into account the different ± cos u and ± sin u factors that multiply the asymptotic expansion of 
the Bessel functions; second, the terms a^kioc) and a2fc+i(a) depend on a, so they change at each level. 

In the sequel, let us fix an integer d > 1, and let us group the different terms of order n~ d l 2 that 
come into play in the sum (|15p : 



If d = 2M is even, then the term corresponding to m = multiplies the Bessel function 
J a (2y/Kz), that we have to expand up to order K~ d , so this term is multiplied by (— l) d cos oj{a). 
The term corresponding to m = 2 will be multiplied by J a+ 2(2\fKz), that we need to expand 
up to order k~ n+1 , because of the {z/Ak) 111 / 2 factor. Taking into account the shift properties of 
the cosine, this is multiplied by (— l)^" 1 coso;(a + 2) = (— 1) cosu)(a). Similarly, all the even 
terms up to m = d will have the same common factor. 

The term corresponding to m = 1 is multiplied by (— 1) sinw (a + 1) = (-l) d - 1 cosa;(a). For 
m = 3, we have (— l) d+1 sincj(a + 3) = (— l) d+1 cosa;(a), and so on. Consequently, we have an 
alternating sum of the form 



2M ,„ / \ 

z \ m/ o I , 1 ^2M-m(a + m) 



S 2M {a,z) = (-l) (j^) P m (a + l,z) m 

m=0 

2M 

= (-4kz)- m ^ {-l) m z m P m {a + 1, z)a 2M - m {a + 



(20) 



for M = 0, 1, . . . 

If d = 2M + 1, then the term corresponding to m = multiplies the Bessel function J a (2y/Kz), 
that we have to expand up to order n~ d , so this term is multiplied by (—1) M+1 sin ui(a). The term 
corresponding to m = 2 will be multiplied by J a+ 2{2y/Kz), that we need to expand up to order 
n~ d+1 , and because of the shift properties of the sine, this is multiplied by — (— 1) A/_1 sinu;(a + 
2) = (— l)^" 1 sinw(a). Similarly, all the even terms up to m = d will have the same common 
factor. 

For m = 1 we have a factor (— 1) M cosuj(a + 1) = (— 1) M smui(a). For m = 3 we obtain 
(_l)M-i cos w ( a + 3) = (— 1) M sino;(a), and so on. 
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Thus, 

2M+1 



S 2M+ l(a,z) = (-l) M+1 ("I) m (^j fln(a + l,z) (2/^M+l-n, 

(21) 

2M+1 v 7 

= {-1) m+1 (4kz)- m - 1 / 2 (-l) m z m P m ( a + l,z)a 2M+l - m (a + m) 

for M = 0, 1, . . . 

Observe that these expansions are in negative powers of k, which is essentially n. A further step 
is needed to retain only the terms up to order n~ d l 2 in each case, but this is simple using symbolic 
computation. 

From (|2(jp and (|21|) . it is possible to compute the coefficients B m (a,z) in (|18|) . Schematically, one 
needs to do the following steps: 

1. Fix the maximum order d. 

2. Generate the polynomials P m (a, z) for m < d, using the recursion (1141) . 

3. Compute the coefficients of the asymptotic expansion of the Bessel functions from (|1T[) up to 
order d. 

4. Compute S m (a, z) for m < d, sum all these terms, expand in inverse powers of n and truncate. 
As an example we will find the first few coefficients. Using Maple, we obtain 

B (a,z) = 1, 



Bi(a, z) 



Az 2 - 12a 2 + 3 
487i ' 
z 3 4a 2 + 11 (4a 2 - l)(4a 2 



Bo(a,z) = 1 z 

y ' ; 288 192 512z 

20a 2 + 187 , 5/2 _ a + 1 3/2 _ (4a 2 - 9)(4a 2 - 25) , 1/2 

31 ' j " 10368 + 23040 48 6144 

(a + l)(4a 2 - 1) (4a 2 - l)(4a 2 - 9)(4a 2 - 25) 

64Z 1 / 2 + 24576z 3 / 2 

„ , , z 6 20a 2 + 391 4 a + 1 , 80a 4 - 584a 2 + 9261 2 

BaIol z) = z H z H z 

y ' ' 497664 829440 576 737280 

(a + l)(4a 2 + 11) z (2a + l)(2a + 3)(2a + 5)(2a + 7)(4a 2 - 32a + 27) 



+ 



384 294912 
(a + l)(4a 2 -l)(4a 2 -9) (4a 2 - l)(4a 2 - 9) (4a 2 - 25) (4a 2 - 49) 
1024z + 1572864^ 2 



Higher order terms become complicated but follow the very same idea. This expansion presents 
two essential advantages with respect to the classical one given by Perron. First, the coefficients are 
still complicated but they can be computed systematically, up to the accuracy desired. Secondly, 
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the original expansion is convergent on the whole complex plane. So, retaining the Bessel functions 
instead of expanding them in negative powers of n, it provides a useful representation of the Laguerre 
polynomials for large degree. 

One difficulty of the previous expansion is that it contains coscj and sinw terms. These terms can 
be grouped together away from [0,oo). Note that 



1 



cosw 



J(2^KZ- 



O 7T 

2 



4 I 



+ e 



-i(2^J~Kz- 



CtTT 

2 4 



- ) 

4 I 



(22) 



If Imz < 0, then the first exponential term dominates. On the other hand, if Imz > 0, then the 
second exponential term dominates. We can write this as follows: 

cosw ~ I e 2v ^ 7 "e ±i (^ L+ ^), ±lmz>0, 

with the understanding that — z = ze Tm if ilmz > 0, i.e. we take the principal argument throughout. 
In a similar way, we can write 



sinw ~ =p 



1 2> /=Si e ±i(^+f) ±i mz>0 , 



2/ 



Now, it is not difficult to check that we can write 

z-/ 2 - 1 / 4 cos W ~ h-zy^-^e 2 ^, 



-a/2-1/4 



' sinw ~ ±—(—z)~ " 



-a/2-1/4 2yf=KZ 



(23) 



ilmz > 0, 



again with the same criterion to understand — z in the complex plane. So we can pull out a common 
factor in the whole expansion: 

Theorem 2 Let a > — 1. The Laguerre polynomial L^\z) admits the following asymptotic expansion 
as n — » oo : 

"d-l 



1 T(n + a + 1) 



2^ 



r «/2-l/4 e V=Si 



Y,B m {a,z)n- m l 2 + 0{n- d l 2 ) 



,m=0 



(24) 



for some coefficients B m (a,z) independent of n. The error term is uniform for z in bounded sets of 
C\ [0, oo) and the coefficients B m (a,z) are related to the original ones B m (a,z) in the following way: 



B2m(a,z) = B 2m (a,z) 
B 2m +i(a, z) = ±iB 2m +i(a, z), ±Imz > 0. 



(25) 



It is important to take into account that this modified expansion is only valid for z in bounded sets 
of C\ [0, oo). On the positive real axis it is not possible to pull out a common factor from the sines and 
cosines, since the argument is real and thus both exponential terms are oscillatory and comparable in 
size. 



Remark 1 When z < we consider arg z = tt, so —z = ze 
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Remark 2 It is possible to remove the ± in the complex plane by writing the odd coefficients in terms 
of the variable —z and using the fact that z m l 2 — (_ z ^m/2 e ±mm/2 ^ ±Imz > and m G Z. For 
instance, 

* , , 4z 2 - 12a 2 + 3 

Bi(q, z) = . , (26) 

T/iis extra step compensates the different sign in i?2 m +i(a, z) in the upper and lower half-plane, 
but one needs to check each power of z within £?2m+i(a, z) individually. 

It is not complicated to recover the standard Perron expansion from the previous one. If we 
compare (|24p with ©, we observe that we just need to expand the terms in the prefactor in inverse 
powers of n and combine them with the coefficients B m (a, z). In this way we have 



■J—z(a _|_ \\ „ 
Ci(a,z) = - K - J -B (a,z)+B 1 (a,z) 



C (a,z) = B (a,z), 

—z(a 4- 
2 

Az 2 - 24(a + l)z - 12a 2 + 3 
~ 48 v /3 i ' 

n( \ (Q+l)(l-2a + z(g + l)) > V=i(g + 1) ■» . 

C 2 (a, z) = B (a, z) H Bi(a, z) + B 2 (a, z) 

o z 

z 3 (a + l)z 2 (20a 2 + 48a + 13)z (2a - l)(2a - 3)(a + 1) (4a 2 - l)(4a 2 - 9) 
~ ~288 + 24 192 32 512 

(27) 

Observe that C±(a,z) coincides with the known coefficient in the Perron expansion, see |18j and 
[19J. The next coefficient C 2 (a,z) agrees with the one given in [U p. 3306], making the changes 
a i — y — a and z t— > — z and reexpanding the prefactor therein, which is written in terms of the shifted 
variable m = n + 1. 



3 Asymptotics of ratios of Laguerre polynomials 

As an application, the above expansion that we have obtained allows us to determine the outer 
asymptotic behavior as n — > oo of arbitrary ratios of Laguerre polynomials with greater accuracy than 
the formulas available in the literature. We write the previous expansion as follows: 

Lt\z) = ft\z) (j2B m (a,z)n- m / 2 + 0(n- d / 2 )) , 



\m=0 



with the prefactor 



f£\*) = ^r T{n + J + ^ e zl \-<n,a)z)-«l 2 -^e 2 ^. (28) 

Note that we emphasize that k depends both on n and a, since we want to consider different degree 
and different parameter of the Laguerre polynomials. Thus 



r(«) ( r \ A") ( \ d-1 
L n {Z) J n {Z) k=Q 
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where 



= , ^ T(n + j + a + l) T{n+l) «(n^)WV 4 o^^jy^^ 
f$P\z) r(n + j + l) T(n + /3 + 1) K (n + j, a)*/ 2 * 1 /* 



Notice that we need different coefficients Bk(a, z) in the numerator, as a result of using degree n+j 
tead of n. However, this is easily computed with the same Maple procedure, setting k = n+j + ^^-. 
The Dfc(a, /3, z) coefficients can be computed in an automatic way. For instance, the first ones are 

D (a,f3,z) = 1, 

f3 2 -a 2 

D 1 {a,P,z) = ±i 4 , ±lmz>0, ^9) 

(f3 2 -q 2 )(3a 2 + 9f3 2 -4z 2 -9) 

D 2 {a,(3,z) = — . 

ybz 

Bearing in mind that z~ 1 / 2 = =pi(— z) _1//2 when ilmz > 0, so we can write 

B 2 -a 2 

D 1 (a,P,z) = ?- 1 =. (30) 
4v— z 



This is the type of modification mentioned in Remark [2] before, rewriting the coefficient in terms 
of —z instead of z, and it becomes more complicated to be implemented for higher terms. 

The prefactor can be expanded in negative powers of n as a convolution of two series. Firstly, the 
ratio of exponential functions is straightforward to compute symbolically. Then 



Lemma 1 We have 

«(n + j,a)«/a+i/* =n 2 1^ A ^^> 

v ' m=0 

where for m > 



^, a ,, )H ^)"f; ( _ ir -,(f t l)(»- + J-t)(^Tj . m 



fc=0 

Proof. We note that 

k(/3 )/V2+1/4 /? + 2j+a + 1 X- 

w 2 1H — ~ — 1 + 



n 2 



fc I \ 2n J ' \ m J V 2n 

k=0 v 7 v 7 m=0 v 7 v 

h V V ( 2 + 'A (^) h (rn-k + Z-\\ k (2J + a + ^ 

2^ h ) \ 2n ) \ m-k { ' \ 2n 



n 2 

m=0 fc=0 

oo 

a 

n 2 

m=0 

where the coefficients A m (j, a, B) are given by ([3"T]) . ■ 

The ratio of Gamma functions can also be expanded as follows: 



(32) 
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Lemma 2 We have 

T(n + j + a + l) j+a _ ^ G m {a,f3,j) 

— — — ~ n J ^ p > , re — > oo, 

T(n + 8 + l) ^ n m 

m=0 

where the coefficients G m (a, (3, j) can be expressed in terms of generalized Bernoulli polynomials: 

G m {a, 0, j) = (j + a -P) B ^-P+i) {j + a + i). 
\ m J 

Proof. See formulas HJ 5.11.13, 5.11.17]. ■ 

Notice that, unfortunately, these generalized Bernoulli polynomials are not directly implemented 
in Maple. However, they are the coefficients of the following generating function: 

00 ±n 

-Y.B^ix)- |t|<27r, (33) 



e t — 1 J ^— ' re! 

' n=0 

see [H 24.16.1]. It is straightforward to compute this expansion symbolically, select the coefficients 
and construct the generalized Bernoulli polynomials from there. 
As a consequence, 



Corollary 1 We have 



r(n + 1) ' n m ' \m/ re r 

m=0 m=0 



r(n + j + l) ,^ G ro (0,0,i) /i\B^ +1) (j + l 



re — )• oo. 



Therefore, this ratio can be computed as a special case of the previous one, setting a = (3 = 0. 

Theorem 3 Let a,/3> —1, and z € C \ [0, oo), i/ie ratio of arbitrary Laguerre polynomials has an 
asymptotic expansion 

r( Q ) ( 7 \ p-a oo 

«~(~) 2 EM«,Aj>K m/2 , (34) 

(2) m=0 

where the first coefficients are 

U (a,P,j,z) = 1, 

/3 2 - a 2 + 2z(/3 -a- 2j) 



(6j 2 +6(a-/3)j + a 2 + 2/3 2 -3a/3)z (/3 2 - a 2 + 2a - l)j ( 35 ) 



U 2 (a,f3,j,z) = — h 

+ 



12 4 

(a 2 - /3 2 - 2a - 2/3 - l)(/3 - a) (a 2 + 3/3 2 - 3)(a 2 - /3 2 ) 



8 32z 
As particular cases, we get the following: 
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If we set /3 = a + 1 and j = 0, we obtain 
and if /3 = a + 2 and j = 0, 



consistently with [9j §3, Lemma 2]. 
If we set /3 = a, then 

C/ (a,a,J» = 1, Ui(a,<x,j,z) = V^zj, U 2 {a,a,j,z) 



so 



*&(*) 1 + vE!- 4 

again in accordance with [9J §3, Lemma 1]. 



a 1 



2 



(2a - 2jz - l)j 



- + 0(n" 3 / 2 ), 
n 



(36) 



(37) 



(38) 
(39) 
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